Wave function statistics and multifractality at the spin quantum Hall transition 
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The statistical properties of wave functions at tfie critical point of tlie spin quantum Hall tran- 
sition are studied. The main emphasis is put onto determination of the spectrum of multifractal 
exponents Ag governing the scaling of moments ~ j^-id-i^q ^j^jj ^jjg gygtem size L and the 

spatial decay of wave function correlations. Two- and three-point correlation functions are calcu- 
lated analytically by means of mapping onto the classical percolation, yielding the values A2 — —1/4 
and A3 = —3/4. The multifractality spectrum obtained from numerical simulations is given with 
a good accuracy by the parabolic approximation A, ~ q{l — q)/8 but shows detectable deviations. 
We also study statistics of the two-point conductance g, in particular, the spectrum of exponents Xq 
characterizing the scaling of the moments {g'')- Relations between the spectra of critical exponents 
of wave functions (A,), conductances (Xq), and Green functions at the localization transition with 
a critical density of states are discussed. 



I. INTRODUCTION 



In the framework of the random matrix theory pio- 
neered by Wigner [Q and Dyson |^ , the statistical prop- 
erties of spectra of complex systems are described by 
random matrix ensembles. Within the Dyson's classifica- 
tion, three symmetry classes are distinguished (orthogo- 
nal, symplectic, and unitary), depending on whether the 
system is invariant under the time-reversal transforma- 
tion and on its spin. It has been understood that this 
classification is very general and applies to a great va- 
riety of physically distinct systems (see for a recent 
review). 

While the Dyson's classification is complete for the 
bulk of the spectrum, more symmetry classes may arise 
in the vicinity of a special point on the energy axis. 
Such non-standard symmetry classes have attracted a 
considerable research attention during the last decade. 
One group of them is formed by three chiral ensembles 
describing the spectrum of a massless Dirac opera- 
tor near zero energy. The same symmetry is shared by 
tight-binding models with purely off-diagonal disorder at 
the band center More recently, four more symme- 
try classes were identified , which characterize a dirty 
superconductor or a mesoscopic superconductor-normal 
metal system. The Hamiltonian matrix has in this case 
an additional block structure in the particle-hole space 
induced by the form of mean-field Bogoliubov-de Gennes 
equations for a superconductor. It was argued |Q that 
the extended classification scheme including 10 classes 
(three Wigner-Dyson, three chiral, and four Bogoliubov- 
de Gennes) is complete. 

The classification of random matrix ensembles can be 
equally well applied to disordered electronic systems. 
In particular, two-dimensional systems of non-standard 
classes are of large interest, in view of their relevance to 
high-Tc superconductors, which have an unconventional 
(d-wave) symmetry of the order parameter and therefore 



possess low-energy quasiparticle excitations. In this pa- 
per, we will consider a system of class C, which is the 
Bogoliubov-de Gennes class with broken time-reversal 
but preserved spin rotation invariance. The correspond- 
ing Hamiltonian satisfies the symmetry H* — —UyHay 
(with Uy the Pauli matrix in the particle-hole space) and 
has the block structure 



H 



h 

A* 



(1) 



where h = h'^ and A = A-^. 

Similarly to the conventional Wigner-Dyson unitary 
class, a two-dimensional system of class C undergoes a 
transition between the phases with different quantized 
values of the Hall conductivity ^-pX]]. More precisely, 
since the quasiparticle charge is not conserved in a su- 
perconductor, one is led to consider the spin conductivity 
determining the spin current as a response to the gradi- 
ent of the Zeeman magnetic field. The quantization of 
the Hall component of the spin conductivity tensor was 
named the spin quantum Hall (SQH) effect. It was shown 
Q that the SQH effect can be realized in superconductors 
with d.j.2_y2 + idxy pairing symmetry explored in recent 
literature 

While the SQH transition shares many common fea- 
tures with its normal counterpart, it is qualitatively dif- 
ferent as concerns the behavior of the density of states 
(DOS) at criticality: while the DOS is uncritical for the 
conventional quantum Hall (QH) transition, it vanishes 
at the SQH critical point. A network model describing 
the SQH transition was constructed in and critical 
exponents for the scaling of the localization length were 
determined numerically. In a mapping onto a super- 
symmetric spin chain was performed, providing an al- 
ternative method for the numerical study of the critical 
behavior. Remarkably, some exact analytical results for 
this problem have been obtained by mapping onto the 
classical percolation [ p^jjl^p^ . Specifically, it was found 
that the DOS scales as p(e) ~ e^/*", while the average 
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product of the retarded and advanced Green functions 
n(r,r') = (Gi?(r,r')GA(r',r)) (referred to as the diffu- 
sion propagator, or the diffuson) and the average two- 
point conductance {g{r, r')) at e = fall off as |r— r'|^^/^. 

It is known that critical wave functions at the con- 
ventional QH transition have multifractal nature ||l^,|l5| . 
Recently, there has been a growth of activity in the direc- 
tion of quantitative characterization of the corresponding 
spectrum of fractal dimensions |16-20|. Zirnbaucr p^ ] 
and Bhaseen et al. |18| proposed a certain supersymmet- 
ric (T-model with a Wess-Zumino-Novikov-Witten term 
(in two slightly different versions) as a candidate for the 
conformal field theory of the QH critical point. The the- 
ory implies an exactly parabolic form of the multifrac- 
tality spectrum. This was confirmed by a thorough nu- 
merical study of the wave function statistics at the QH 
transition ||20[| . 

The aim of this paper is to study the wave function 
statistics at the SQH critical point. We will demonstrate 
that the exponents A2 and A3 governing the scaling of 
the second and third moments of the wave function in- 
tensity (see Sec. || for the formal definition) can be calcu- 
lated exactly by analytical means. Quite surprisingly, we 
find that the index rj = — A2 characterizing the spatial 
decay of the wave function correlations is equal to 1/4, 
in contrast to the r~^/^ decay of the diffusion propaga- 
tor. This leads us to a general analysis of relations be- 
tween different critical exponents characterizing the wave 
function statistics in the qualitatively new situation of 
the localization transition with a critical DOS. We com- 
plement our analytical results by numerical simulations, 
which allow us, in particular, to investigate whether the 
multifractality spectrum of the SQH critical point is ex- 
actly parabolic or not. The answer to this question, as 
well as the exact values of A2 and A3 we have found, is of 
central importance for identification of conformal theory 
of the SQH transition, which is the issue of a considerable 
research interest at present [^- ^. Some of our results 
were reported in a brief form in ||24| . 

The article is organized as follows. In Sec. |l| we re- 
mind the reader of some basic concepts related to the 
multifractality of critical wave functions. In Sec. ^ we 
describe the network model of class C and use it to calcu- 
late numerically the DOS at the critical point of the SQH 
transition. In Sec. ^ we present an analytical calculation 
which involves a mapping onto the percolation theory and 
allows us to calculate the averages of products of two and 
three Green's functions and thus the exponents A2 and 
A3. Section ^ is devoted to a numerical evaluation of 
the full multifractal spectrum A^. This allows us not 

but also to 



[V 



only to check the analytical results of Sec. 
investigate whether the spectrum is exactly parabolic (as 
for the conventional QH critical point) or not. In Sec. ^ 
we present a numerical study of statistical properties of 
the two-point conductance. We further include a scaling 
analysis of the relation between the multifractal spectra 



of the two-point conductance and of the wave functions 
at a critical point with a vanishing DOS. These analytical 
arguments clarify the connection between the numerical 
findings of Sec. VI and the results of Sec. IV, ^ on the 
wave function multifractality. Finally, Sec. VII contains 



a summary of our results and a brief discussion of some 
remaining open problems. 



II. WAVE FUNCTION MULTIFRACTALITY IN 
SYSTEMS WITH NON-CRITICAL DOS 

Multifractality of wave functions ip{r) is known to be 
a hallmark of the localization transition. It has been ex- 
tensively studied in the context of conventional Ander- 
son and quantum Hall (QH) transitions with non- critical 
DOS (see ||r^ , [l5| , ^ and references therein), and we re- 
mind the reader of some basic results. Multifractality is 
characterized by a set of exponents 



Tg = d{q -1) + Aq 



(2) 



(d is the spatial dimensionality) describing the scaling of 
the moments of |?A^(r)| with the system size L, 



(|V(r)p^)^L 



-d—T„ 



(3) 



Anomalous dimensions Ag distinguish a critical point 
from the metallic phase and determine the scale de- 
pendence of wave function correlations. Among them, 
A2 = —rj plays the most prominent role, governing the 
spatial correlations of the "intensity" , 



L^^{\i;\v)i?{v')\)^{\v-v'\/L)- 



(4) 



Equation (^ can be obtained from (||) by using the fact 
that the wave function amplitudes become essentially un- 
correlatedat |r— r'| ^ L. Scaling behavior of higher order 
spatial correlations, (ri)V'^''^ (r2) • • -^^^"(rn)!) can 

be found in a similar way. Correlations of two differ- 
ent (but close in energy) eigenfunctions and the diffusion 
propagator n(r,r';c.) - (Gi+^(r, r')G^(r', r)) {G^-^ 
are retarded and advanced Green functions) possess the 
same scaling properties. 



L 



^'^dV'KOVKr')!) 

L2<i(^,(r)V;*(r)V'*(r')V,(r')) 
2H(r, 



P 



3 

,r'; w) 



(5) 



where cj = — e^-, L^^ ^ (pcj)^^/'', p is the density of 
states, and |r — r'| < L^. In two dimensions the multi- 
fractal spectrum Aq plays a key role in the identification 
of the conformal field theory governing the critical point, 
which led to growing interest in the eigenfunction statis- 
tics at the QH transition [l^[2C|]. 

Applying naively these results to the SQH transition, 
one would conclude that the r~^/^ scaling of the diffu- 
sion propagator found in ||l^ implies rj = 1/2. However, 
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we show below that this conclusion is incorrect. This 
demonstrates that one should be cautious when trying to 
apply the relations between critical exponents obtained 
for systems with a non-critical DOS to those with a crit- 
ical one (like the SQH transition), as will be discussed in 
Sec. and Sec. ^, 



III. 



NETWORK MODEL AND THE DENSITY OF 
STATES 



As a model of the SQH system, we use the SU(2) ver- 
sion Q of the Chalker-Coddington network describing 
the QH transition |^^. Dynamics of the wave function 
defined on edges of the network is governed by a unitary 
evolution operator U. At each node of the network the 
scattering from two incoming into two outgoing links is 
described by a matrix 



S = 



cos tf 
— sin( 



sm t 
cos ( 



(6) 



with — tt/A corresponding to the critical point. Each 
realization of the network is characterized by a set of 
random 2x2 spin matrices Ue associated with all edges 
e of the network. In view of (|l|), U satisfies the sym- 
metry U = (TyWcTy, implying that Ue G SU(2). Diago- 
nalizing U for a square network of the size L x L yields 



4L eigenfunctions ipia{e) and eigenvalues e * 
i = 1, 2, . . . , and a = 1, 2 is the spin index. 



where 




FIG. 1. Scaling plot of the density of states for system 
sizes L = 16(<>), 32(n), 96(o). Dashed and dotted lines indi- 



es), 5 = l/27rL^/* 



cate power laws (dashed: e ' , dotted: 
denotes the level spacing at e = 0. Inset: same data on a 
linear scale and the result from the random matrix theory 
(solid curve). 

We begin by displaying in Fig. ^ the numerically cal- 
culated DOS p(e) for different system sizes L. It is seen 



that after a proper rescaling all data collapse onto a sin- 
gle curve. Specifically, the energy axis is rescaled to 
t/5, where 5 oc i"^/** is the level spacing at e = 0. 
(This scaling of 5 is related to the critical behavior 
of DOS p{e) ^ e^/'' discussed below via the condition 
p{5)5 ^ Xjl? ) The scale invariance of pit) at critical- 
ity is reminiscent of the analogous property of the level 
statistics at the conventional Anderson or QH transition 
(see for a review). At e ^ 5 the critical DOS scales 
as pit) ~ e^^^, in agreement with analytical predictions 
pO| . On the other hand, at e ~ (5 one observes an oscil- 
latory structure qualitatively analogous to the behavior 
found in the random matrix theory (RMT) for the class 
C §. 

Let us note that, strictly speaking, deviations of DOS 
from the RMT at e i5 are not parametrically small. 
On the other hand, the numerically found DOS follows 
very closely the RMT curve for two oscillation periods. 
In other words, the energy scale below which the RMT 
works (the effective Thouless energy), while being para- 
metrically of order (5, turns out to be several times larger. 
This indicates that there is a numerical smallness in the 
problem, and the SQH critical point shows "close-to- 
metal" features (similar to the Anderson transition in 
2 -|- e dimensions with small e). The small value 1/7 of 
the DOS exponent is another manifestation of the same 
fact. 

The states with energies t^ b are localized with the 
localization length t~'^l'^ ||l^. For smallest energies 
t ^ b the correlation length is of the order of the sys- 
tem size. In view of their critical nature, these states are 
expected to be multifractal, L^''(|7/'iQ(e)p'') ~ L"^". For 
e 3> (5 the multifractality holds within a region of the ex- 
tent (outside which the wave function is exponentially 
small); hence 



-2(g-l)-A, ^ 



(7) 



By the same token, spatial correlations are expected to 
be governed by the multifractality on scales below . In 
particular, we have for correlations of two different eigen- 
functions with energies e^, tj t 



iMV'.a(e)V,73(e')P>-(^/?e) 



A, 



< 



6 



(8) 



(r is the distance between e and e'), and similarly for 
higher-order correlators. In Sec. [V and ^ we will demon- 
strate the multifractality explicitly and calculate the ex- 
ponents Aq. 



IV. TWO-AND THREE-POINT CORRELATION 
FUNCTIONS: MAPPING ONTO PERCOLATION 
PROBLEM 

In this section, we present an analytical calculation 
of two-point and three-point correlation functions, which 
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allows us to find the fractal dimensions A2 and A3. We 
use the mapping onto the classical percolation, follow- 
ing the approach of Q , and demonstrate that it can be 
extended on products of two and three Green's functions. 



A. Two-point functions 

Consider a correlation function of two wavefunctions, 
P(e',e;ei,e2) = VL(e)Vja(e)V»/3(e')V'*/3(e') 

ijaf3 

X 5{ei-e,)5{t2-ej)), (9) 

where e, e' are two different edges of the network. Intro- 
ducing the Green function 

G(e',e;z) = {e'\{l - zU)-^\e) 

(which is a 2 X 2 matrix in the spin space), we express 
(i) as 

P(e', e; ei, 63) - (27r)-2(Tr[Gfl(e', e; e''^) - GA{e' , e; e'^^ )] 
x[Gfl(e,e';e"^)-GA(e,e';e^^=)]), (10) 

where Gr^a are retarded and advanced Green functions, 
Gfl,A(e',e;e"i) = G(e', e; e*(<^i±'°)). We will calculate 
( |l0| ) at zero energy, £1^2 — > 0, but finite level broadening, 
±iO zti'j with 7^1. The scaling behavior of the cor- 
relation function (^ at £1,62 ~ e can then be obtained 
by substituting e for 7. We thus need to calculate 

V{e\ e; 7) = (2^)-2(Tr[G(e', e; z) - G(e', e; z'^)] 

X [G(e,e';z)-G(e,e';z-i)]), (11) 

with a real z = e"''' < 1. By the same token, in order to 
understand the scaling properties of another correlator 
of two wave functions, 

P(e',e;ei,e2) = (5]|^^a(e)n^,/3(e')P 

ijaf3 

X S{ei - e,;)(5(e2 - Cj)), 



(12) 



we will consider the correlation function 



D{e\ e; 7) = (2^)-2(Tr[G(e, e; z) - G(e, e; z~')] 



X Tr[G(e', e'; z) - G(e', e'; z-^)]). (13) 



As discussed in the end of Sec. Ill, the scaling behavior of 
( pi] ) and ( |l3| ) at r <C (where r is the distance between e 
and e') is governed by the multifractal properties of wave 
functions (specifically, by the exponent A2). The general 
strategy of calculation of the correlation functions (11), 
( p3| ) is analogous to that used in |l^] for the one-point 
function TrG(e, e;2;). Therefore, we outline only briefly 
those steps which generalize directly the calculation in 
[fi3[, and concentrate on qualitatively new aspects. 



The Green functions in ( [III ) , ( [T^ are straightforwardly 
represented in the form of a sum over paths 



G{e,e';z)= ^ 

paths e'- 



(14) 



where Sj is the corresponding matrix element (cos 9, sin 6, 
or — sinO) of the 5'-matrix between the edges Cj and e^+i. 
Equation ( |l^ ) generates a convergent expansion in pow- 
ers of z when \z\ < 1; otherwise the identity 



Gt(e,e';z) = l-,5ee' -G(e',e; (z*)'^) 



(15) 



is to be used (in all our calculations z is real, so that 
z* — z). As shown below, each of the double sums over 
paths obtained by substituting ([l|), (|l|) in (|ll|) or (|l|) 
can be reduced to a single sum over classical paths (hulls) 
in the percolation problem. This remarkable reduction 
crucially relies on the following two statements: 

1. Only paths visiting each edge of the network either 
or 2 times are to be taken into account; contribu- 
tions of all the remaining paths sum up to zero, 

2. Using the statement 1, it is easy to see that each 
node may be visited 0, 2, or 4 times. The second 
statement concerns the nodes visited four times. 
As illustrated in Fig. ||, there are three possibilities 
how this may happen; the corresponding contribu- 
tions have weights (i) cos"* 9, (ii) sin"* 9, and (iii) 
— sin^ 9 cos^ 9 from the scattering matrix at this 
node. The statement is that one can equivalently 
take into account only the contributions (i) and (ii) 
with the weights cos^ 9 and sin^ 9, respectively. 



i) 



ii) 



iii) 

s c 



-s c 



FIG. 2. Possible configurations of paths passing four times 
through a network node. The symbols c and ±s denote the 
elements cos'^ 6, ± sin^ of the S-matrix at the node. 



In Ref. both statements were proven for the case 
of the average of a single Green function (G(e, e; z)). We 
show below that they are valid for all the two-point func- 
tions entering (^l|), (p^), as well as for averaged products 
of three Green's functions (considered in Sec. IV B). Let 
us emphasize that such a generalization is far from triv- 
ial. This point is well illustrated by the fact that prod- 
ucts of four (or more) Green functions determining the 
exponents A, with g = 4, 5, . . . can not be mapped onto 
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the percolation within our approach (see Sec. IV C| and 
Appendix) . 

We now proceed by proving the statement 1 . It is con- 
venient for us to recall first the corresponding proof for 
the case of a single Green function, (TrG(e, e; z)), consid- 
ered in For an arbitrary edge / the paths entering 
( p^ ) can be classified according to the number k of times 
they pass through /. The contribution of paths with 
k ^ has the form 



Y,{T^-B[UfAifJ)f 



(16) 



fc=i 



where Bi is a sum over all paths going from f to e and 
then from e to /, and A{ f, /) denotes a sum over paths 
which begin and end on / and do not return to / in 
between. Since A{f, /) is a linear combination of SU(2) 
matrices with real coefficients, it can be represented as 
Aif,.f) = |A(/,/)|i(/,/), where A{fJ) e 5[/(2) and 
\A{f,f)\ is a real number. After a change of the inte- 
gration variable, UfA{f, /) — ^ Uf, Eq. ( [l^ ) then reduces 
to 



(17) 



fe=i 



Since SU{2) matrices can be represented a.s U = 
exp(ian<T), with a real a and a unit vector n {ai are 
the Pauli matrices), one finds 



U — cos ka ■ 1 + i sin ka ncr. 



(18) 



The SU(2) invariant measure is {2/n) da sin^ a f dn, 
where dn is the conventional measure on the sphere. 
Therefore, for an integer k 



{U'') = Ck ■ 1, 



Cfc 



1, 

1 

2 ' 



k = 
/c = 2, 



(19) 



0, otherwise. 



Substituting (|l9|) in (|T^), one finds that only the term 
with k — 2 survives, which completes the proof of the 
the statement 1 for the case of an average of one Green 
function. 

We turn now to the products of two Green functions. 
Consider 

{Ti-G{e,e';z)G{e",e"';z)) (20) 

(we will need below both cases e" = e', e'" = e, and 
e" = e, e'" — e'). Using (p^, we classify the contribu- 
tions to (EG) according to the numbers of returns fci, k2 



to the edge / for the corresponding two paths. We want 
to show that only the contributions with fci -t- ^2 = 0, 2 
are to be taken into account. If one of ki is zero, the 
proof is obtained in the same way as for a single Green 
function (see above). We thus consider the remaining 
contributions, which are of the following form: 



oc 

E 

fci ,k2 = l 



(Tri?i [UfAif, f)f^ B2 [UfAif, f)f^ ) , (21) 



where Bi is a sum over the paths f e and e'" —>■ /, 
and i?2 is a sum over the paths / e" and e' — > /. Per- 
forming the variable change UfA{f,f) — > Uf as before, 
we get 



^ {TTB,U^'B2U^')\A{fJt'+''\ (22) 

fci,fc2 = l 

Using ([l^), we calculate now the average over Uf in (|2^): 
{TyBiU^j^ B2Uf) = TrSiB2(cosfciacosfc2a)a 

Tr _Bi(Tii?20'i(sinfciQ;sinfc2a)a 



= iTr(SiS2 + ^ E BiaiB2(Ti)ck^+k2 

i 

+ ^Tr{BiB2 - ^^Bia,B2a.)cfc,_fe, 



biCki 



+k2 



■k2- 



(23) 



The only property of the factors 6i, &2 which is impor- 
tant for us at this stage is that they are independent of 
ki, k2. The sum (p^ is therefore reduced to the form 



J2 ib^Ck,+k2+b2Ck,-k2)\Aifjt^+'^. (24) 

While the first term in brackets is non-zero only for 
ki + k2 = 2 {i.e. fci = ^2 = 1) as required, the second 
one seems to spoil the proof. Let us perform, however, a 
summation over ki at fixed ki + k2 = k. Using Eq. ([is|), 
we find then that the coefficients in the second term can- 
cel for any even k > 4 (for odd k all terms are trivially 
zero): 

Cfci„fe2 = Cfc-2 + Cfc_4 + . . . + C_(i._2) 



ki+k2 = k 



= C2 + Co + C_2 = 0. 



(25) 



Therefore, only the term with ki — k2 — 1 survives in 
the sum (p5), which completes the proof. ^ 



^The correlation function (TrG(e, e; 2:)TrG(e', e'; z)) is ana- 
lyzed in the same way, yielding again a sum of the type (p^, 
so that our argument remains valid. 
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Applying now the statement 2, the proof of which is 
given in Appendix, we represent each node as a super- 
position of contributions of the types (i) and (ii) (Fig. ||) 
with weights cos^ 6 and sin^ 9, equal to 1/2 at the SQH 
critical point. The network is then reduced to a weighted 
sum over all its possible decompositions in a set of closed 
loops (such that each edge belongs to exactly one loop). 
These loops can be viewed ||l^,|l3| as hulls of the bond 
percolation problem. Non-zero contributions to the cor- 
relation function ( pO| ) come from pairs of paths retracing 
exactly twice a loop or a part of it. This yields for z < 1 

(TrG(e', e; z)G(e, e'; z)) = (TrG(e', e; z-^)G{e, e'; z'^)) 
= -2^F(e',e;7V)z2^, (26) 

N 

(TrG2(e',e;z)) = (TrG2(e, e'; z"!)) 

= -^Pi(e',e;7V)z2^, (3?) 

N 

where P{e',e;N) and Pi{e',e;N) are probabilities that 
the edges e and e' belong to the same loop of the length 
(resp. with the length A^ of the part corresponding 
to the motion from e to e'). Furthermore, to calculate 
the correlation function (TrG(e', e; z)G(e, e'; z"-^)) enter- 
ing ([1^), we apply the identity (|l5| ) to the second Green 
function and then use the property 

(TrG(e', e; z)G\e' , e; z)) = -2{Ti-G^{e', e; z)) (28) 

following from the SU(2) symmetry. As a result, we find 



(TrG(e',e;z)G(e,e';z"i)) = -2 ^ Fi(e', e; A^)^ 



2JV 



N 



(29) 



and, combining ( |26D and (|29|), 

TT^Die', e; 7) = i Y^^Piie', e; N) + P, (e, e'; A^)]; 



.2N 



N 



J2P{e',e;N)^ 



2N 



(30) 



N 



Equations (|26D, (p9[), ( pO[) express the quantum corre- 
lation functions entering (pi) in terms of purely classical 
quantities P(e', e; A^) and Pi(e', e; A^). To analyze the re- 
sults, we recall some facts from the percolation theory. 
It is known that the fractal dimension of the percolation 
hulls is 7/4 implying (see ||2^ for a recent discussion) 
that P and Pi scale as 

P(e', e, N), Pi(e', e, N) ~ N-^'^-^'^ , r < A^^/v 

(31) 

and fall off exponentially fast at r ^ N^^'^ , where r is 
the distance between e and e'. This yields for the corre- 
lation functions in (|2^) and ( p9| ) (which we abbreviate as 
{GrGr), (GaGa), (GrGa)) 



{GrGr) = (GaGa) 0^ (GrGa) - 



7 



^4/7 



(32) 



in full agreement with the scaling argument of . How- 
ever, these leading order terms cancel in (B3) since 



P{e', e, TV) = ^ Pi(e', e, N) = P(e', e), (33) 



N 



N 



where P(e',e) is the probability that the edges e and e' 
belong to the same loop. The result is non-zero due to 
the factors z"^^ only, implying that relevant A^ are now 
A^ - 7^1, so that {{Gr -Ga){Gr-Ga)) scales differ- 
ently compared to (|32|), 

P(e', e; 7) = ^ ^) " ^)Kl " e"'"^^) 

N 

^P(r, 7-^)7-' -(C7O"'/', ^<C7- (34) 
Using now the definition (||) of V and the DOS scaling, 



P(e) 



.1/7 



we find for r < 



L4(^*Je)V',.(e)^.Me')V'*Me')) ^ ('^/fj 



-1/4 



(35) 



The correlation function (^3|) is calculated in a similar 
way. The results for the {GrGr), (GaGa), and (GrGa) 
terms in (n3) have the form 



(TrG(e, e; z)TrG(e', e'; z)) 

= 4-2^[P(e;Ar)-f P(e';A^)]; 



2^ 



N 



NN 



N 



^ P_ (e, e'; A^, iV')z2(^+A^') + ^ ^),2Ar^ (gg^ 

(37) 



(TrG(e,e;z^i)TrG(e',e';z-i)) 
= 5] P_ (e, e'; TV, 7V')z^(^+^') -|- ^ P(e, e'; A^) 



,2Ar 



WW 



.27V 



(TrG(e, e; z)TrG(e', e'; z^^)) - 2 ^ P(e'; Af);: 

AT 

- 5] P_ (e, e'; A^, iV')z2(^'+^^') - ^ P(e, e'; Ar)z2^, (38) 



N 



where P(e; A^) is the probability that e belongs to a loop 
of the length A^, while P-(e, e'; A^, A^') is the probabihty 
that e and e' belong to different loops of the length A^ 
and A^', respectively. A larger number of terms in (|3^)- 
(^) as compared to (^), ( ^9|) is because of two reasons. 
First, there is a unit matrix contribution of a "path of 
zero length" to the expansion (|lj) of the Green function 
G(e,e; z). Second, e and e' may now belong to different 
loops and still give a finite contribution, since each of the 
two paths will retrace twice the corresponding loop. 



Combining (36)-(p3) and using the identities 
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P{e,e';N)+J2P-ie,e;N,N')=P{e;N), 

N' 

N 

we get for the correlation function jl^ ) 

^2p(e',e;7) = 5]P(e,e';7V)(l-z2A^) 

N 



(39) 
(40) 



(41) 



AT TV' 



We see again that a.t z = e'' ~ 1 the rcsuh is zero, and 
that at smaU 7 it is dominated by N ^ 7^^. Using 
Eq. (|l|) and 

P_ (e, e'; TV, N') - P(e; iV)P(e'; iV') - Ar-8/^(Ar')"*^^ 

(42) 



we find that the first term in (^) is ~ ^. 



-1/4-1/4 



at 



— 1/2 

7" ^ while the second one is ~ and thus can 

be neglected. Therefore, we find that "D shows the same 
scaling behavior as V [see Eq. (^], 



P(e',e;7)^5]P(e,e';7V)(l 



-2N-f\ 



N 

,^-l/V-l/4^ 



(43) 



In other words, the wave function correlator 
{\'^iai^)'^jl3(.^')ll with ei,ej ~ e scales at r < in the 
same way as (|35|), 

i'(IV^L(e)V'|/3(e')l> - p-^(e)P(e',e;7 e) ^ (r/e.)-^/^ 

(44) 



Both Eqs. (35) and (44) imply that the fractal exponent 

1 



?7 = -A2 



(45) 



at variance with what one might naively expect from 
the r~^^'^ scaling of the diffusion propagator (GuGa), 
Eq. (H). 



B. Three-point functions 

We consider now averaged products of three Green 
functions, analogous to the two-point functions ( [III ) and 

V{e, e', e"; 7) - (2^)-3(Tr[G(e, e'; z) ~ G(e, e'; z'^)] 
X [G(e',e";z) -G(e',e";z-i)] 
X [G(e", e; z) - G(e", e; z^^)]), (46) 

P(e, e', e"; 7) = (2^)-3(Tr[G(e, e; z) - G(e, e; z"!)] 
X Tr[G(e',e';z) -G(e',e';z-i)] 
X Tr[G(e", e"; z) - G(e", e"; z-^)]}. (47) 



The key role in the calculation of ( [46[ ) and (47) is played 
by the proofs of applicability of the statements 1 and 2 



(Sec. [V A) to the products of three Green functions. De- 
tails of these proofs are given in Appendix. After the two 
statements are applied and the network is reduced to a 



sum over its loop decompositions (as in Sec. IV A), the 
correlation functions are calculated straightforwardly. In 
particular, we find for the averaged products of three 
Green functions entering (pq) 



(TrG(e, e'; z)G(e', e"; z)G(e", e; z)} 

= - ^[3P(e, e', e"; N) + P(e", e', e; N)\: 



-2N 



N 



(TrG(e, e'; z-i)G(e', e"; z-i)G(e", e; z"!)} 
= e', e"; iV) + 3P(e", e', e; A^)].; 



-2N 



N 



(TrG(e, e'; z)G(e', e"; z)G(e", e; z^^)) 
--2^Fi(e,e',e";Ar); 



-2^ 



N 



(TrG(e, e'; z-i)G(e', e"; z-i)G(e", e; z)) 
= 2^Pi(e",e',e;Af)^ 



-2Af 



(48) 
(49) 
(50) 
(51) 



where P(e, e', e"; A^) is the probability that the edges e, 
e' and e" belong to the same loop of the length A^, with 
e' lying on the path from e" to e, while Pi{e,e',e";N) 
is the same probability but with A^ being the length of 
the segment from e" to e. Combining Eqs. (^8|)-(pTl), 
we express the correlation function (^6|) in terms of the 
classical probabilities P and Pi. Remarkably, the situ- 
ation is qualitatively different as compared to the cal- 



culation of two-point functions (Sec. IV A): the leading 
terms in (^)-(^l|) do not cancel in the expression for 
2?(e, e', e"; 7). We can thus simply set 7 = (z = 1), 
which yields 

(27r)32?(e, e', e"; 7) ^ 2[P(e, e', e") + P(e", e', e)], 

r<^ij, (52) 

where P(e, e', e") = P(e, e', e"; A^) is the probability 
for e, e', and e" to belong to the same loop with the ori- 
entation e ^ e' ^ e" ^ e, and r is the characteristic 
scale of the distances between e, e', and e". 

The correlation function ( |47| ) is calculated in the same 
way, and the results are qualitatively similar. We thus 
skip intermediate formulas and only present the final re- 
sult, 

(27r)3p(e, e', e"; 7) c 8[P(e, e', e") + Pie", e', e)], 

r « C^,, (53) 

which differs from (p2h by an overall factor of 4 only. 
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Using any of the equations (p2|), (p3|), we can deter- 
mine the fractal exponent A3. In analogy with (|3l|), the 
probability for the edges e, e', and e" separated by dis- 
tances ~ r to belong to the same loop (percolation hull) 
of a length N scales as 



P(e, e', e"; N) ~ iV^*^/ V^/^, r < N^" 



(54) 



and is exponentially small for r ^ N^/"^ . Summing over 
TV, we thus get 

P(e,e',e") ~r-3/4. (55) 



Substituting this in Eqs. (|52|), ( |53| ) and expressing P and 
T) in terms of wave functions in analogy with the two- 
point functions (ph, (p2|), we find for r < 



L6(V>,„(e)V'*^(e')V',Me')V';^(e")^fc^(e")V'L(e)) 



i'(|^»(e)V,73(e')V'fc7(e")r) 



J.-3/4 



p3(e) 



-3/4 



Therefore, the exponent A3 is equal to 

3 

A3 = — . 

^ 4 



(56) 



(57) 



C. Discussion 

The situation we encountered while calculating two- 
and three-point functions is qualitatively different from 
what happens at conventional localization transitions. 
Specifically, in the conventional case average products 
of only retarded or only advanced Green functions are 
negligible compared to mixed averages containing both 
Gr and Ga. e.g. {GrGr), (GaGa) < (GrGa). For 
this reason, the wave function correlators, which are pro- 
portional to {{Gr — Ga){Gr — Ga)}, are determined by 
{GrGa) (and similarly for higher moments). In contrast, 
we have found in the SQH case that the correlators of 
the {GrGr) (or {GaGa)) type are approximately equal 
to {GrGa) and cancel it in the leading order (so that 
{{Gr — Ga){Gr — Ga)) scales differently). Evaluation of 
three-point functions made the overall picture even more 
complex: while we obtained again an identical scaling of, 
say, {GrGrGr) and {GrGrGa) correlators, this time 
the cancelation was not complete, and the correlation 
function {{Gr, — Ga){Gr — Ga){Gr — Ga)) scaled in the 
same way. 



To shed more light on the reason for these different 
types of scaling behavior, it is instructive to reverse the 
logic and to examine how the diffuson scaling ds^) can 
be obtained from wave function correlations (^5[). It is 
straightforward to express the zero-energy diffusion prop- 
agator in terms of the correlation function I?(e', e; ei, €2) 
defined in Eq. (|l^), 



n(e', 



(1 



(TrG«(e',e;l)G^(e,e';l)) 

—^P(e,e; 61,62) 



)(1 



(58) 



As discussed in Sec. [V A , I?(e', e; 61, 62) scales with the 
distance r — \e' — e| and the energy 61^2 ~ 6 as follows 



P(e',e;6i,62)~(r/C.)^^e 



(59) 



where A2 — —1/4, Xp 
of DOS defined by p{e) 



-4/7 



1/4 is the scaling dimension 
and = e-i/(2-^p) ^ 
(For r T^{e' , e; ei, 62) is exponentially small.) 
Substituting ( p9| ) in (|5^), we see that if 2xp -f A2 > 0, 
which is the case for the SQH transition, the energy in- 
tegral in ( ^8|) is dominated by 61,2 ~ e(?'), where e(r) is 
defined by ^^(j.) r {i.e. e{r) ~ r^(^-^p) = r"''/'^), and 
can be estimated as 



n(e',e)^P(e',e;ei,62)|e,, 



(60) 



in full agreement with the exact result (|34).^ This is in 
a stark contrast with the case of a conventional local- 
ization (Anderson or QH) transition, when the diffusion 
propagator 11 (or any other correlation function of the 
{GrGa) type) depends in a singular way on the infrared 
cutoff set by L^^ , see the last line of Eq. (^) . On the other 
hand, Eq. ( |60| ) has a familiar form of a two-point corre- 
lator in a conformal field theory (or, more generally, in 
field-theoretical description of standard critical phenom- 
ena), where {Oi{ri)Oi{r2)) scales as |ri — r2\~^^\ with 
Xi being the scaling dimension of the operator Oi. 

Generalization to higher moments is straightforward. 
We define a wave function correlation function 

p(«)(ei,...,e,;6i,...,6,) = (27r)-« 

X {Tt[{Gr - Ga) (ei , 62 ; e^^ ) {Gr - G^) (e2 , 63; e^'^ ) 

X...X (Gfl-GA)(e„ei;e'^')]), (61) 

and a set of {G . . .G) correlation functions, 

nif...s,(ei,...,e5;6i,...,6,) 

= (TrG,, (ei, 62; • ■ • G,,(e„ ei; e^")), (62) 



^ Since the integral (^^ is determined by the upper cutoff 
e(r) (and not by the vicinity of e = 0), this calculation ap- 
plies not only to 11 = {GrGa), but equally well to {GrGr) 
and {GaGa), in agreement with (B2). 
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where Sj = R or A. Assuming that all distances between 
the points Ci are ~ r, we have in analogy with (p9), 



(63) 



Writing for ni'-'..s, a spectral representation of the type 
( p8| ) , we see that the integrals are determined by the up- 
per limit e ~ e(r) = provided 



r{q) = qxp 



Aq > 0. 



(9) 



Under this condition, wc find that 11, 
pendent of the indices Si and scales as 



(64) 



is in fact inde- 



n(9) 



-qXp 



(65) 



For larger q, when qXp + < 0, the energy integrals are 
dominated by the vicinity of e ~ 0. Consequently, the 
correlation functions nii^...,s, start to depend in a singu- 
lar way on the infrared cutoff' (^e) and are expected to 
scale in the same way as 'D^'^\ Eq. ( |63| ) (with a numer- 
ical prefactor depending on indices Si), similarly to the 
conventional Anderson localization transition. 

The value of q separating the two regimes is thus de- 
termined by the equation qXp + Ag = 0. For the SQH 
transition {xp = 1/4) its solution is, in view of Eq. (57), 
g = 3. Remarkably, this is also the largest value of q for 
which the mapping onto percolation described above still 
works (see Appendix) . We believe that this is not a mere 
coincidence. Indeed, within this mapping average prod- 



ucts n. 



(9) 



of Green functions are expressed in terms of 
probabilities of the percolation theory, and are therefore 
of order unity for r ~ 1. On the other hand, Eq. (^) 
yields, in the regime qXp + Aq < 0, a result which is 
much larger than unity at r 1, ^ 1 and diverges 
in the absence of the infrared cutoff, — > oo. We see 
no way how such a behavior might be produced by the 
percolation theory. 

Finally, we discuss a relation between our consideration 
and the field-theoretical approach to the wave-function 
multifractality ||2^ |3^,|l^,^. In the rcnormalization- 
group language, T{q) defined by Eq. (|6j) are scal- 
ing dimensions of operators of the type O'-''^ ~ 
ipsi'ipli ■ ■ ■'4'sq'4'l' J where V'j V"^ ^'^^ electronic fields. Av- 
eraged products of Green functions are expressed as cor- 
relation functions of the corresponding operators in 
particular, (B2|) takes the form 



(Tr0W^(e2)OW3(e3)...OW,(ei)}. (66) 



To calculate the scaling behavior of such correlation func- 
tions, one applies the operator product expansion (OPE) 
p0|"p2[ . Generically, the identity operator will be among 
those generated by the OPE. Moreover, under the con- 
dition r{q) > [Eq. p3)] it will be the most relevant 



operator and will dominate the expansion, leading to the 
gap scaling 11^''' 7--9r(i)^ agreement with (^5|). On 
the other hand, if T{q) < 0, the operator O^"^^ will give a 
dominant contribution to OPE, leading to a multifractal 
type of scaling, n(«) oc r^^^^i) (r/^,)i'(9) , as in Eq. (|6|). 
What is, however, non-trivial from this point of view, 
is that the scaling of the wave function correlator ( |6l| ) 
has the multifractal form (|6^) independently of the sign 
of r{q). This means that in the regime T{q) > the 
leading (gap scaling) terms ( pq ) cancel in the particular 
combination of the functions IT'^-' corresponding to I?'''' , 
and subleading terms determine the result (|6^). A simi- 
lar cancellation of leading scaling terms in the context of 
classical percolation was recently discussed in [bsl . 



V. WAVE FUNCTION STATISTICS: 
NUMERICAL RESULTS 

A. Multifractality spectrum 

The analytical treatment of Sec. |^ yielded results for 
the anomalous dimensions at two distinct values of q, 
A2 — —1/4 and A3 = —3/4. In order to obtain more 
complete information about the wave function statistics, 
namely the multifractality spectrum at arbitrary q, we 
have performed numerical simulations. A question we are 
particularly interested in is whether or not the spectrum 
is exactly parabolic. A definite answer on this question 
will imply, along with exact values of A2 and A3, an im- 
portant constraint on the conformal theory of the critical 
point, which is a subject of current research [pT]-p3|. 

Before we come to the presentation of our findings, 
we give a few remarks about technical aspects of our 
numerics. We compute wave functions by numerically 
diagonalizing the AL^ x 4L^ unitary time evolution op- 
erator U of the Chalker-Coddington network described 
in Sec. 



Ill 



Using advanced sparse matrix packages [p3[ , 
we selectively calculate only states with energies in the 
vicinity of e = 0, which are critical over the whole ex- 
tent of the system (^^ ~ L). Specifically, we consider, for 
each realization of the network, four lowest eigenstates 
{i.e. with eigenvalues e"'*^ closest to unity). The num- 
ber of wave functions in a statistical ensemble we obtain 
this way ranges from about 10^ for i = 16 to 2 • 10"' for 
L = 384. 

To determine the multifractality spectrum r^, we cal- 
culate for each wave function ipi the generalized inverse 
participation ratio (IPR) 



2q 



(67) 



and analyze the scaling of the average (Pq) with the sys- 
tem size L. The data can be fitted very well by the power 
law 
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(P,)=c,(2L)--'. (68) 



of q we studied, 



To demonstrate this, we show in Fig. |^ the system size 
dependence of {Pq){2LY'' , with Tq obtained from the fit. 
The plot is organized in such a way that a pure power 
law (^8|) would correspond to a horizontal line. This 
kind of plot is very sensitive to any corrections to a pure 
power-law behavior of {Pq). Since no systematic cur- 
vature is observed, corrections to scaling are extremely 
small. This allows us to determine the anomalous dimen- 
sions Aq Tg + 2(1 — q) with great accuracy. 
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FIG. 3. Scaling of the average IPR with the system size L 
for several values oi q = 0.5(o), 1.5(0), 2(a), 2.5(<i), 3(V), 
3.5(l>), 4(n). The system size dependence of the amplitude 
Cq{L) = {Pq)(2LY'' is presented, with = 2{q - 1) + Aq 



shown in Fig. ^ The scattering of the data is due to the lim- 
ited size of the statistical ensemble used. The solid line is a 
guide to the eye corresponding to the vanishing of finite size 
corrections {cq{L) = const). 



are shown by a solid 
0. We choose to plot 



The obtained results for Aq 
line in the upper panel of Fig. 
Aq/g(l — q), since this would give a constant for an ex- 
actly parabolic spectrum, which is uniquely determined 
by 77, Aq = r]q{l — q)/2. According to our analytical cal- 



culations (Sec. IV), Aq/q{\ — q) is equal to 1/8 for both 



q ~ 2 and q ~ 3; this value is marked by the dashed line 
in the figure. It is seen that the numerical results agree 
perfectly well with the analytical findings at g = 2 and 
q — 3. Furthermore, the parabolic dependence may serve 
as a numerically good approximation in the whole range 



(69) 



Nevertheless, we believe that Eq. (^) is not exact. In- 
deed, at < g < 2 the numerically found show clear 
deviations from exact parabolicity (69), which are of the 
order of 10% near q = 0. Since this is precisely the regime 
in which finite-size effects have been found to be very 
weak and Aq was determined with a high accuracy, we 
interpret the observed deviations as very strong evidence 
for nonparabolicity of the exact multifractal spectrum of 
the SQH transition. In particular, the deviation of the 
limiting value Aq/q{l - q)\q^o = 0.137± 0.003 from 1/8 
well exceeds the estimated numerical uncertainty. 



0.14 




FIG. 4. Upper panel: Anomalous dimension A, (solid line) 
describing the scaling of the average IPR {Pq)- The functional 
form Ag/qjl — q) highlights the deviation from exact parabol- 
icity (p9|) indicated by the dashed line. The circles correspond 
to the exponent A, obtained from the scaling of the typical 
value P,'yP. 

Lower panel: Singularity spectrum f(a). numerical results 
(solid line) and the parabolic approximation ( |70| ) (dashed 
line) are shown. The inset depicts a magnification of the 
apex region; the deviations from (^o|) correspond to the en- 
hancement of Aq/q{l — q) near q — in the upper graph. 



ptyp 
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We also calculated typical inverse participation ratios, 
= exp(lnPq) and the corresponding dimensions 
2{q — 1) + Aq.^ It follows from the general analysis 



of the wave function multifractality |34| that f„ 



Tq for 



^In our earlier publications p4| , po|j36[ we used the symbol Tq 
to characterize the scaling of the typical value Pq^^, and fq 
for the average (Pq). In the present paper we have chosen to 
interchange the notations. 
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Q < 9c, where qc corresponds to the zero a_ of the singu- 
larity spectrum /(a) (defined below). In the present case 
we find from the data qc = 3.9 ± 0.1 (the parabolic 
approximation (^9|) would imply qc = 4). For q > Qc 
the average (Pg) is determined by rare realizations, and 
Tg < fq. Furthermore, already for q smaller than but 
close to qc, finite-size corrections to P^^p become large 
pO| , leading to large errors in determination of Tq. For 
the SQH problem, we find that the scaling of P*^p ex- 
hibit small finite-size corrections as long q < 2.5, so that 
the corresponding exponents fq can be found with a high 
accuracy. The results are shown by circles in Fig. ^ (up- 
per panel) and are in full agreement with the values of 
Tq obtained from the scaling of (Pq). For larger q {q > 3) 
the finite-size corrections to P^^p (which we estimate to 
be ^ L~y with y x 0.4) become appreciable, strongly 
reducing the accuracy of determination of fq. 

The lower panel of Fig. || depicts the singularity spec- 
trum f{a) obtained by a numerical Legendre transform 
of the scaling dimension Tq, f{aq) = qaq — Tq with 
ctq — dTq/dq. The dashed line represents the parabohc 



(70) 



const X q^{q — 1)^ 



(71) 



approximation corresponding to (p9| 



/(a) - 2 - 



[a - a^f 



"0-2 = 1/ 



4(ao - 2) ' 

We see again that the parabolic approximation is numer- 
ically rather good; nevertheless, it is not exact. Devia- 
tions from ( [70| ) are demonstrated in the inset which shows 
an enlarged view of a region around the maximum of 
/(a). The deviation of ao - 2 = 0.137 ± 0.003 from 1/8 
corresponds to non-parabolicity of Tq discussed above. 

B. IPR fluctuations 

We devote the remainder of this section to a brief dis- 
cussion of the IPR distribution function V{Pq), specif- 
ically, its evolution with the system size L and depen- 
dence on q. In analogy with Anderson and quantum 
Hall transitions studied earlier |34|j2C| , |35| -|37|| , we expect 
the distribution V{Pq) to become scale-invariant in the 
large-L limit. Figure o demonstrates that this is indeed 
the case. It represents the evolution of the distribution 
of lnP2 with the system size L. The mean of the distri- 
bution is shifted as —f2 In L. Apart from small statistical 
fluctuations at the largest system sizes, a clear tendency 
towards an asymptotic form is observed. To characterize 
the width of the distribution P(lnPg), we calculate the 
variance cr^ — var(lnPg), as shown in the inset of Fig. |^ 
for g = 2. The results extrapolated to L — > oo (the finite- 
size corrections are again of the type L'^ with y « 0.4) 
are presented in Fig. ^ The behavior of cr^ is qualita- 
tively similar to that found for other localization transi- 
tions. A somewhat unusual feature of the SQH transition 
is that in a rather broad range < g < 3 the variance ct^ 
is remarkably well described by the formula 



which has been derived for a metallic system ||3^,|25|], or 
for the Anderson transition with a weak multifractality, 
e.g. in 2 + e dimensions ||3^ , |36|] . In the latter case this 
formula is valid for q <^ qc. The accuracy of Eq. ( [7l| ) is 
one more manifestation of the "close-to-metal" character 
of the SQH critical point already mentioned in Sec. IH, 
which leads to a relatively large value of (?c — 4. At 
larger q, the behavior of aq becomes linear (as was also 
found for the conventional Anderson transition ||3^ , ^7|| ), 
in agreement with the theoretical prediction Uq ~ q/ qc for 
q ^ qc [Q. This is because in this regime the distribu- 
tion V(Pq) is dominated by a slowly decaying power-law 



tail V{Pq)(X Pq 



where Xq — qdq for g > ge p4 



Q. 



2 




FIG. 5. Distribution function 'P(lnP2) for system sizes 
L = 16, 24, 32, 48, 64, 96, 128, 192, 256, 384 (from right to 
left). For values of L large enough, the form of the 
distribution becomes independent of L. Inset: Width 
a2(L) = ((lnP2 - (lnP2))^)^/^ of the distribution P(lnP2) 
versus the system size L. Some scattering of the data for the 
largest system sizes is due to the limited number of samples. 
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(75) 



FIG. 6. R.m.s. deviation a, = {{InPq - (lnPq))Y^^ of 
the IPR logarithm. The dashed line is a fit to Eq. (|7l|), the 
dotted Une corresponds to the asymptotic limit a, = q/qc 
with qc = 3.9. The inset shows an enlarged view of the low-g 
region. 



VI. STATISTICS OF TWO-POINT 
CONDUCTANCES 

So far, we have investigated properties of an open sys- 
tem. To define the two-terminal conductance g, one 
opens the system by attaching two leads. According to 
the Landauer-Biittiker formula, g = TitH, where t is the 
transmission matrix between the leads. In the framework 
of network models, the transmission matrix determining 
the two-point conductance between the edges e and e' 
has the form ||T^ 



t = {e'\{l-UP)-^U\e), 



(72) 



where P = 1 — |e)(e| — |e')(e'| projects out the states on 
the edges e and e'. 

Statistics of the two-point conductance g{Y,Y') has 
been extensively studied, both analytically and numer- 
ically, for the conventional QH transition, exemplifying 
a localization transition with non-critical DOS. It was 
shown that the moments {g'^{r',r)) obey a power-law 
scaling [|9|jll0, 

(5nr,r'))'-|r-rr^', (73) 



with a set of exponents Xq related to Aq |19| , P0| , 



2Ai/2 , 



q<l/2 
(7 > 1/2. 



(74) 



For the SQH critical point, only the average conductance 
{g{e' , e)) has been considered previously. Gruzberg, Lud- 
wig, and Read pO| found that 



Beamond, Cardy, and Chalker ||T^ used the mapping 
onto percolation to calculate {g{e', e)) at the band center 
e = (z = 1), with the result 



(g(e',e)) = 2P(e',e), 



(76) 



with P(e',e) as defined in Eq. (|33|). Comparing 
( [76| ) with (p9|), we see that {g{e',e)) is equal (up 
to the sign) to the diffusion propagator n(e',e) = 
(TrGfl(e',e;l)G^(e,e';l)). 

In this section, we will study statistical properties of 
g{e',e) at the SQH transition. Note that though the 
definition of {g{e' , e)) reminds closely that of the diffu- 
sion propagator 11 (e', e), the identical scaling of the both 
quantities is not at all self-evident. In contrast, they scale 
differently at conventional localization transitions, as can 
be easily seen by comparing Eqs. ([73|), ( [zi] ) with (^. It 
is worthwhile to remind the reader the physical reason 
for this difference (see also a related discussion in p6[). 
The product (G'^(r', r)G';^(r, r')) has a meaning of the 
particle density (or, in an optical analogy, the radiation 
intensity) at a point r' induced by a source inserted into 
the system at a point r. In an infinite system at critical- 
ity this quantity turns out to be infrared divergent: if a 
source is switched on at a time t = 0, the detected inten- 
sity will increase with time without saturation, since the 
radiation cannot propagate away fast enough. Therefore, 
in order to make (G^G'^) finite, one needs to allow the 
propagating wave to get out of the system, i.e. to intro- 
duce absorption. One possibility is to make the absorp- 
tion weak but uniform over the whole system , leading 
to (Gf+,^(r',r)G^_,^(r,r')) =n(r',r; 2^7), which is the 
same as introducing a small uniform level broadening 7 
(or equivalently, a small frequency lu with an analytical 
continuation to the imaginary axis, lu — 2ij). Alterna- 
tively, one can allow for a particle to be absorbed at the 
points r and r' only, but with a probability of order unity, 
yielding the two-point conductance g{r',r). Clearly, two 
definitions are essentially different (which is already obvi- 
ous from the very fact that 11 depends on 7, diverging in 
the limit 7^0, while g does not require any parameter 
like 7 and is bounded, g < I). Therefore, the different 
scaling behavior of 11, Eq. (^ and (g), Eq. (fz^), is not 
surprising. 

Returning to the SQH transition, we are thus naturally 
led to a question: why do H and (g) scale identically in 
this case? The reason is that the zero-energy diffusion 
propagator H(e',e) is in fact defined at 7 = {i.e there 
is no need to introduce absorption or a finite frequency 
to regularize it), see Eqs. ( ^ and (^), (60). This can 
be traced back to vanishing of DOS at e = 0. It is not 
a surprise that in this situation, when the absorption is 
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irrelevant, n(e',e) and {g{e',e)) (which only differ in the 
way the absorption is incorporated) scale in the same 
way. 

Let us consider now higher moments 

([TrG(e',e;e-^)G(e,e';e-'r). 

Applying the consideration of Sec. IV C| , we find that the 
absorption (7) remains irrelevant provided 



2qxp + A25 > 0, 



(77) 



with the result 
([TrG(e',e;, 



-T)G(e, e'; e^)]«) - r-^ii^^+^''\ (78) 



For the SQH case, the condition (|7^) implies q < 3/2. 
(We make an assumption that our consideration, which 
is strictly speaking performed for integer q, remains valid 
for intermediate, non-integer values of q.) According to 
the above argument, ((7*) scales in this regime in the same 
way (|78|), so that (see also [p3[). 



Xq = 2qxp + 2Ag, 
and, using Xp = 1/4 and Eq. (|6S|), 



Xg ^ q{3 - q)/A. 



(79) 



(80) 



Note that, in contrast to TrGG, the two-terminal con- 
ductance g is bounded from above, g < 2 (the factor two 
is due to spin summation and is not essential). Phys- 
ically, it simply means that for such rare realizations 
when TrGG is large, g is limited by the contact resis- 
tance. It follows that the exponent for {g'^) should be a 
non-decreasing function of q. In other words, the expo- 
nent Xq saturates after reaching its maximum at some 
go- We find from ( |80| ) go — 3/2; for larger q the exponent 
saturates at the value Xg^gg ~ Xgg ~ 9/16 (these mo- 
ments are determined by the probability to find g 1). 

Equation ( |80| ) implies, in particular, a normal distribu- 
tion of In g at r ^ 1 with the average {\ng{r)) = —Xt Inr 
and the variance var[ln(7(r)] = 61nr, where Xt ~ 3/4 
and b ~ 1/2. These values correspond to the parabolic 
approximation (|69|); more accurate predictions can be 
obtained by using the nmnerical results for Ag, 



Xt=XQ = 2xp 
b 



2A'^ 



2^\q — '^iX p 
0.58 



2(ao - 2) 



0.774, (81) 
(82) 



(here a prime denotes the derivative with respect to q). 

We turn now a to a numerical study of the two-point 
conductance. While we did not attempt a high-precision 
numerical determination of the spectrum of correspond- 
ing exponents Xg (as presented in Sec. VA for the mul- 
tifractal spectrum of wave functions), we have verified 
some of the key predictions of the above analytical con- 
siderations. Figure ^illustrates evolution of the distribu- 
tion function V{g) with the distance r between the con- 
tacts; it is seen that at sufficiently large r the distribution 



becomes log-normal as expected. In Fig. g we show the 
scaling of the average {g) and the typical gtyp = exp(ln (7) 
values of the two-point conductance, along with anal- 
ogous quantities (|G|^) and |G|typ — exp(ln|Gp) for a 
closed system, |Gp = -TrG(e', e; l)G(e, e'; 1). For the 
average values, {g) and (|Gp), the numerics fully confirm 
the theoretical results ([76| ), ( p9| ) telling us that the both 
quantities scale as r~^/^iid, moreover, are equal to each 
other. A non-trivial character of the equality (g) = (|Gp) 
is well illustrated by the data for typical quantities: gtyp 
and |G|jyp are not equal. Nevertheless, they are found to 
share a common scaling: gtyp, |G|jyp ~ r^^*, confirming 
our arguments presented above. Furthermore, the nu- 
merically obtained value of the exponent, Xt ~ 3/4, is in 
agreement with the theoretical prediction (|8l|). 
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FIG. 7. Distribution of the two point conductance at 3 dif- 
ferent distances between the contacts, r — 5.7(*), 14.1(-|-), 
133(v); the system size is L = 196. The dashed line indi- 
cates a log-normal fit with parameters {Ing} = —4.72 and 
var(lngi) = 3.15. 
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FIG. 8. Scaling of the two-point conductance with distance 
r between the contacts: average value (empty symbols), (g), 
and typical value (filled symbols), gtyp — exp(ln g), in systems 
of sizes L — 128(a) and L — 196(o). Also shown is scaling of 
the two-point Green function, (IGj^) and |G|typ = exp(ln |Gp) 
(L = 128(A), L = 196(o)). The lines correspond to the r"^/^ 
(dotted) and r""^/* (dashed) power laws. Deviations from 
power-law scaling at large values of r are due to the finite 
system size. 



by the conformal field theory methods. Note that our 
results are against the proposal of Ref. [^, where the 
result Aq = q{l — q)/A was obtained. Apparently, this 
indicates that the theory considered in and obtained 
p2[ from a particular network model with fine-tuned cou- 
plings, does not belong to the SQH universality class. 
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VII. SUMMARY 

Let us summarize the main results of the paper. 

1. We have extended the mapping of the SQH net- 
work model onto the classical percolation and cal- 
culated two- and three-point correlation functions 
at the SQH transition. This allowed us to deter- 
mine analytically the fractal exponents A2 and A3 
governing the scaling of the second and third mo- 
ments of the wave function intensity, with the re- 
sults A2 = -?7 = -1/4 and A3 = -3/4. 

2. We have performed a thorough numerical study of 
the multifractal spectrum A^. The obtained spec- 
trum is given with a good accuracy by the parabolic 
law ( |69| ) but shows clear deviations from parabol- 
icity. Fig. |. 

3. Statistical properties of generalized inverse partic- 
ipation ratios Pq at the SQH transition are sim- 
ilar to those found earlier for other localization 
transitions. In particular, the distribution function 
V{Pq) becomes scale-invariant in the limit of large 
system size. 



4. We have analyzed statistics of the two-point con- 
ductance g at the localization transition with a crit- 
ical density of states. Specifically, we have pre- 
sented scaling arguments which link the exponents 
Xq governing the spatial decay of (5') to the wave- 
function multifractality spectrum Ag, see Eq. (79). 
This yields, in particular, for the typical conduc- 



-Xt 



with 



tance at the SQH critical point gtyp 
Xt ~ 3/4 (see Eq. ( |8l| ) for a more accurate value), 
as confirmed by numerical simulations. 

In recent years, a considerable progress has been 
made in understanding of conformal field theories re- 
lated to problems of two-dimensional fermions subject 
to quenched disorder 
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1J,|21| g3|,|0|-|4^. In partic- 
ular, a relation between the wave function multifractal- 
ity in two-dimensional disordered systems and the opera- 
tor content of corresponding conformal field theories has 
been discussed in a number of publications |3^,^,|l^,^ . 
It remains an open question whether the multifractal ex- 
ponents Aq, Xq for the SQH transition can be computed 



Discussions and correspondence with J.T. Chalker, 
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and D.G. Polyakov are gratefully acknowledged. This 
work was supported by the Schwerpunktprogramm 
" Quanten-Hall-Systeme" and the SFB 195 der Deutschen 
Forschungsgemeinschaft . 



APPENDIX A: PROOFS TO THE MAPPING 
ONTO PERCOLATION 



In Sec. IV A two statements were formulated which al- 



low us to calculate averaged products of two (Sec. [V A) 
and three (Sec. IV B) Green functions. Here we give some 
more details on the proofs of these statements. In the end 
of the Appendix we will also explain why our calculation 
cannot be extended on products of g > 4 Green functions 
(and thus on higher moments of wave functions). 



Statement 1 

The first statement says that only paths visiting each 
node or 2 times are to be considered. Its proof for the 
case of two-point functions was given in Sec. [V A| . The 
analysis of the case of three-point functions (considered 
in Sec. IV B) goes along similar lines, and we present 
its brief outline only. In analogy with (^2|), we have to 
consider an expression of the type 



E 



(TrSi C/^i ^2 C/^" B3 C/^^ ) I A (/, /) I ''^1 



(Al) 



where ki is the number of returns of the i-th path 
{i — 1,2,3) to the edge/. Performing averaging over 
Uf as in (p3|), we cast (Al) into the following form 



00 



/ , iblCki + k2 + k3 + b2Cki+k2-ka + hCki+ka-ki 

^1,^2:^3—1 



+biCk2+ka-k,)\A{fJt^+''^+''^. 



(A2) 



The first term in curly brackets is trivially zero in view 
of (p^. To demonstrate that remaining terms give 
zero as well, we perform a summation over ki at fixed 
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k = ki + k2 + k^. Indeed, it is not difficult to show by a 
straightforward arithmetics that for an arbitrary k 



Cki+k2-k3 

^1,2,3 — 1,2,...; ki+k2+ks — k 



0. 



(A3) 



Therefore, the sum (A2) is equal to zero, which com- 
pletes the proof of the statement 1 for the three-point 
Green functions. 



Statement 2 

The second statement allows us to reduce the nodes 
visited 4 times to a superposition of contributions (i) and 
(ii) of Fig. ||, with the factors cos^ 9 and sin^ 9, respec- 
tively. We will give the proof for the (most non-trivial) 
case of a product of g = 3 Green functions; the proof for 
q — I and 2 is obtained in the same way. More specifi- 
cally, we will consider the correlation function (^7|); the 
correlator (|4^) is treated analogously. 

Each of three Green functions in (^^ generates a sum 
over closed loops (e e, e' ^ e', and e" — > e", re- 
spectively). For a given lattice node, let us label the 
corresponding incoming edges as (1,2) and the outgoing 
ones as (3,4). We are considering a contribution of paths 
visiting this node in total four times. This generates four 
path segments starting each on one of the edges (3,4) and 
ending on one of the edges (1,2), and not passing through 
any of these edges. We are going to show that for any con- 
figuration of these four segments the statement 2 holds. 
It is easy to see that there exist two essentially different 
types of such configurations (shown in Fig. all others 
can be obtained by permutations of e, e', and e", and/or 
by lattice symmetry operations. 



A'(2,3) 




b) 



A(l,3) 



A(2,3) 



A'(l,4) 




A(l,4) 



A(2,4) 



FIG. 9. Configurations of paths for the correlation function 
( p7[ ) for a node visited four times. 

Consider first the configuration (a) of Fig. |^. We have 
to sum over all ways to connect the four path segments 
by various configurations of scattering events at this node 
shown in Fig. || Specifically, only such connections are 
allowed which generate exactly 3 closed loops, each con- 
taining one of the edges e, e', and e". There are three 



possibilities how this can be done, since the segment that 
does not contain any of the edges e, e', e" can be con- 
nected in a loop with any of the remaining three. In one 
of these cases the configuration of paths at the considered 
node is of the type (ii) of Fig. ||, in two other cases it is 
of the type (iii). We thus get the following contributions: 

(TrC/2A(2,3)[/3 TrC/2A'(2,3)C/3 

xTTUiA{l,4:)U4,UiA'{l,4)Ui) x s'*, (A4) 
(TrC/2A(2, 3)U3 TtUiA{1, A)Ua 

xTrC/2A'(2,3)[/3C/iA'(l,4)t/4) x (-c^s^), (A5) 
(Tr[/2^'(2, 3)C/3 Tri7iA(l, A)Ui 

xTr[/2A(2,3)C73C/iA'(l,4)C/4) x (-c^s^). (A6) 

Here A{2, 3) is a sum over all paths from 3 to 2 pass- 
ing through e, A' (2, 3) is a sum over paths 3 — >• e' — >■ 2, 
A{1, 4) is a sum over paths 4 — > e" 1, and A'{1, 4) is a 
sum over paths 4^1 (Fig. |a). Also, we have denoted 
s = sin 9 and c — cos 9. 

To perform the integration over Ui, we use the follow- 
ing formulas of integration over SU(2) matrices: 



{T:tUViUV2)u = -iTrFiV2, 
{TrUViU^V2)u - ^Trl/i T1V2. 



(A7) 
(A8) 
(A9) 



Here matrices Vi,2 are assumed to be of the form Vi — 
\Vi\Vi, where Vi G SU{2) and \Vi\ is a real number (we 



remind that A(i, j) are exactly of this type, see Sec. IV A ). 

Applying repeatedly the rules (A7)-(|AS|), we 
perform integration over all matrices Ui {i — 
1,...,4) in Eqs. ( |A4| ) -(|A6|) . We find then that all 
three contributions (A4)-(A6) are proportional to 

with coefficients 
The total 



contributions 
Tr^t(2,3)A'(2,3) TrA^XT, 4)^1, 4), 



— ifi^, — ic^s^, and — ic^s^, respectively, 
coefficient is therefore 



I4 I22 ^22 I2 



(AlO) 



We see that the same result would be obtained if we 
would assign the weight to the first contribution 
(which is of the type (ii) of Fig. |^) and discard the re- 
maining two terms (which are of the type (iii)). This 
establishes the validity of the statement 2 with respect 
to the configuration (a) of Fig. |£ . 

The configuration (b) of Fig. 9| is analyzed along the 
same lines. We have again three contributions, one of the 
type (ii) of Fig. || and two of the type (iii) , 



{TrU2A{2,3)U3 TrC/iA(l, 4)C/4 

xTTUiA{l,3)U3U2A{2,A)Ui) x 
(TrC/2 A(2, 3)U3 TtUiA{1, 3)U3 



(All) 
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xTrt/iA(l,4)t/4C/2^(2,4)(74) x (-c^s^), (A12) 
{TtUiA{1, 3)C/3 TtUiA{1, 4)Ui 

xTrU2A{2,3)U3U2A{2,A)Ui) x (-c^s^). (A13) 

After integration over Ui according to the rules 



(|A7|)-(|A9D they all produce an identical structure, 
TrAt(2,3)A(2,4)At(l,4)A(l,3), with the coefficients 
js"^, ^c^s^, and ^c^s^, respectively. Again, retaining 



only the (ii)-type contribution (All) and assigning the 



weight s to it, we would obtain the same result. This 
completes the proof of statement 2 for the three-point 
correlation function (p7|). 



What about g > 3? 

A natural question is whether the present approach 
can be generalized to higher-order correlations of wave 
functions governed by multifractal exponents Aq with 
q > 3. The answer is negative. In fact, both statements 
1 and 2 do not apply (or, in a more careful formulation, 
our proofs fail) for q > 4, as we are going to explain 
in brief. Concerning the statement 1, consider a gener- 
alization of the expression (Al) to q = 4, and choose 
fc = fci + fc2 + fca + fc4 = 4. Obviously, there is just one 
such term (all ki = 1) in the sum, and it is easy to see 
that it is generically non-zero. Therefore, no cancella- 
tion of terms with k > 2 happens in this case, i.e. the 
statement 1 does not work. Turning to the statement 
2, consider e.g. a correlation function £'(e, e', e", e'"; 7) 
analogous to ( ^ ) but containing a product of four traces 
of Green functions. Trying to prove the statement 2, we 
will then have to consider the path configurations very 
similar to those shown in Fig. |9| but with all four paths 
containing one of the edges e, e', e", or e'". At the next 
step the paths should be connected via the scattering 
processes at the node - this time to generate 4 closed 
loops. However, for each of the configurations shown in 
Fig. ^ there is only one way to do this, so that only one 
con tribu ti on wi ll arise in place of three terms ( A4)-([A^) 
or ( All )-( A13 ). Clearly, the statement 2 is not valid in 
this situation. Therefore, the mapping onto the classical 
percolation is not applicable for higher moments, q > 3. 
This is in correspondence with the fact that q ~ 3 sep- 
arates two regimes of qualitatively different behavior of 
correlation functions, as discussed in Sec. IV C. 
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